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A short review of Schrodinger hamiltonians for which the spectral problem has been 
related in the literature to the distribution of the prime numbers is presented here. We 
notice a possible connection between prime numbers and centrifugal inversions in black 
holes and suggest that this remarkable link could be directly studied within trapped 
Bose-Einstein condensates. In addition, when referring to the factorizing operators of 
Pitkanen and Castro and collaborators, we perform a mathematical extension allowing 
a more standard supersymmetric approach. 

PACS number(s): 98.80.Hw, ll.30.Pb 
1. Introduction 

The problem of the nontrivial zeros of Riemann's zeta function, i.e., whether they 
all lie on the line z ~ i + it in the complex plane or not, is a famous unsolved 
mathematical problem in which physics, especially quantum mechanics and chaos 
theory, could have a substantial and rewarding contribution in view of its direct 
connection with the distribution of prime numbers. "'^ The phase of Riemann's zeta 
function could be considered the main mathematical object able to trace this im- 
portant arithmetic distribution. 

In this work, my first goal is to provide a short survey of the hamiltonians 
that so far have been proposed to give hints for a spectral solution (also known as 
the Hilbert-Polya conjecture) of the location of the zeta zeros on the critical line. 
Along the review the reader can encounter several relevant novel implications in this 
highly interesting topic, especially in sections 2 and 4. Since the logderivative of the 
phase shift of the Coulomb repulsive potential with reversed centrifugal barrier is 
a function that approximates very well the Riemann phase for the negative orbital 
number I = —1/4 and because centrifugal inversions have been conjectured in black 
holes a promising connection between prime numbers and black hole solutions in 
general relativity is suggested. Moreover, this could be under direct experimental 
investigation in the not quite far future if sonic black hole configurations could be 
produced in trapped Bose-Einstein condensates. In addition, when addressing the 
factorizing operators of Pitkanen and Castro and collaborators, a simple mathe- 
matical extension of their work is dealt that brings their approach closer to the 
standard supersymmetric methods of quantum mechanics. 
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2. Bhaduri, Khare, Low (1995) ^ 

BKL showed that the density of zeros of Riemann's C function is determined by its 
phase 9{t) 

r fi + ii) 

exp[2i9{t)] = exp{-it\nn)Rr{t) , Rr{t) = ^ _| ' (1) 

^ U 2 / 

where the convention 6{0) = tt has been introduced. Although this phase is smooth, 
i.e., it does not include the jumps by n due to the zeros of the modulus of C, it counts 
well the zeros on the critical line. It can also be expressed as the logderivative of 
the I = ~ \ phase shift (argF (| + ^)) with respect to the distorted Coulomb wave 
of the repulsive Coulomb potential 

dy^ y 

where Z = -i, e = -fp, and k=^. 

The potential in Eq. (2) is also equivalent to an inverted harmonic half-oscillator. 

Of course, such a negative fractional phase shift does not look physical at first glance 
as it points to an attractive centrifugal barrier (i.e., well) but this is precisely the 
case for its application to prime numbers. Even more, I hail this as a very appealing 
feature because of a paper of Cirone et al,^ where the —1/4 correction is called the 
quantum anti- centrifugal force being a metric effect related to the radial derivatives 
in the Laplacian. Furthermore it is well established that a 'centrifugal force reversal' 
occurs for r < 3m in a Schwarzschild spacetime.^ The connection between prime 
numbers and black holes is of course a distingiiished one and of much transcendence. 

It is not hopeless to succeed in touching experimentally the promising con- 
nections alluded above. Especially suitable is the BEC physics. According to 
Garay, sonic black holes could be created in BECs trapped by means of sufficiently 
tight ring-shaped external potentials produced with state-of-the-art or planned 
technology.^ In this case, the sonic black holes are stationary solutions of the Gross- 
Pitaevskii nonlinear equation displaying a boundary from which only a tiny thermal 
quantum sound radiation can escape. Of course, a detailed study of the instablil- 
ity of these sonic horizons is needed as recently emphasized by Leonhardt and 
collaborators.^ 

But most promising, in 2001, Ott et alj reported the production of a BE con- 
densate in a microstructurcd magnetic surface trap. This is remarkable because the 
booming realm of surface-patterning technological processes opens up for the beau- 
tiful physical phase-shift approach to prime numbers. Thus, by means of micro- 
fabricated electrical circuits, precision measurement of scattering properties be- 
comes feasible and interferometric measurements should allow soon a direct access 
of the phase shifts that was never available in the past.^ 



(2) 
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Fig. 1 The BKL Coulomb potential with —1/4 centrifugal correction and the other 

parameters scaled to unity. 

3. Berry and Keating (1999) ^ 

This is in the class of Berry and Keating contributions, who were stimulated by 
ideas and remarks of Connes. According to Berry and Keating, the hamiltonian 



= —ih ( ; " 



Hbk = -ih \^xdx + -j (3) 

is the simplest hermitian operator connected to Riemann's (, though it is only a 
canonically rotated form of the upturned harmonic oscillator — .t^. The simplicity 
results from avoiding the complications of the parabolic cylinder eigenfunctions. 
The linear independent wave functions in the x representation are of the type 

, , ^, A9{±x) 

where is the Heaviside step function. One can also write them in terms of the 
smoothed counting function of the Riemann zeros 

^.(x:E) = °-pH;(;:r:^g-'i , (5) 

where Berry and Keating define A^sm(-B) = 0R{E)/Tr + 1, and 

^fl(i^) = -f logTr + JmlogrQ + iiE^ . (6) 

In addition, Aneva,^° considered chaos quantization conditions for Hbk and their 
geometrical (group) interpretation in an effort to get a discrete spectrum for Hbk- 
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4. Castro et al (2001, 2002) 

Castro and collaborators, based on works of Pitkanen,^^ introduced the simple and 
appealing nonrelativistic supersymmetric quantum mechanical (SUSYQM) schemes 
for the Riemann zeros problem. They also elaborated on p-adic oscillators for the 
same problem. Their first order "ladder" operators are 

d dV{t) , ^ ^ d dVillt) , 

where kc is a constant parameter and V{t) is given by the logarithm of the Gauss- 
Jacobi theta series G(i'), e2^(*) = G(f'), I = const 

OO 

G{t^)= e-^"'*' (8) 

n= — oo 

and V{l/t) 7^ V{t). The resulting pair of susy partner hamiltonians have the same 

real spectrum s(l — s) on both the real line (trivial zeros) and Riemann's critical 
line (nontrivial zeros). The 'coherent' eigenfunctions are written ^s{t) and \l/s(l/t), 
respectively, where 

exp{—t) 



1 — exp(— t) 



(9) 



The Mellin transform of the Gauss-Jacobi theta series is the zeta function ({s) and 
the inner product of two eigenfunctions, say ^'^i and ^'s2, can be written in terms 
of (^{a{sl +S2) + b— |), where a and b can be expressed in terms of ki, k2 and I. The 
zeta zeros are in one-to-one correspondence with the 'orthogonality' conditions. 

Here I will show that there is a standard way of doing SUSYQM for this impor- 
tant case. Using a relationship of Gauss that was mentioned by Castro 

G{l/x) = x^'^G{x) , (10) 

one can write the above ladder operators in the form 

d dV(t) , , ^ d dV(t) , 

where k2 ~ ki + j. 
This leads to 

Ha^a = D^D^^A = i-d^t + {ki - k2)dt + dtV + + (fci + k2)V + /ei/ea]*^ (12) 
and 

Hb^b = DiD2^B = [-d^t + (ki - k2)dt - dtV + V^ + (fci + k2)V + /ci/c2]*b , (13) 
where dt = -j^&ndV =^ lnG{t^). 
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Using now the gauge transformation "9 a.b = '^A,Bt''^^ one can eliminate the first 
derivative dt and get the Schroedinger equations 

dttpA,B + Sa,B-4'A,B = ^kiM'^A^B , (14) 

where ^^kxM — (^^^j^)^ ~ ^1^2 plays the role of eigenvalue and Sa,b are in the 
position of Schroedinger susy partner potentials 

Sa,b = ±dtV + + {ki + k2)V . (15) 



5. Chadan and Musette (1993) 

These authors used 



(p ['(^ + 1) - /cm(a;;5,a,/5) 
where the centrifugal correction of Chadan and Musette is 



(16) 



f / p\ —9 S /3 

_ Jcm[x;g,a,p) _ l^{x) l{x) l{x)ll(x) , 

^ cm 9 — 9 • v-^'/ 

X X 

In Eq. (17) L = log(l/a;), LL = loglog(l/a;), ; = 1, a = -27/4, /3 = -3/2 and g is 
a coupling constant. 

Chadan and Musette proposed the above rather complicated singular potential in 
a closed interval [0,i?] and Dirichlet boundary conditions at both ends. They gave 
arguments that the spectrum in the coupling constant 5 = i + ^ (t is the imaginary 
part on the critical line), which is real and discrete, with g„ > 1/4, coincides 
approximately with the nontrivial Riemann zeros when R — e^^'^/'^. 

We note that this is a so-called Sturmian quantum problem, i.e., a quantization 
problem in the coupling constant of the potential. A very detailed analysis of this 
singular hamiltonian and a generalization thereof from the point of view of inverse 
scattering and s-wave Jost functions has been performed in the important work of 
Khuri.^'^ 




Fig. 2 The centrifugal correction Fcm of Chadan and Musette for R < 0.015 and g < 50. 
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6. Joffily (2003): zeros of s-wave Jost functions 

The first to relate the nontrivial zeta zeros to the complex poles of a scattering 
matrix were Pavlov and Fadeev,^^ for the case of a particle on a surface of negative 
curvature. 

Joffily, on the other hand, suggested a 'potential scattering' Hilbert-Polya con- 
jecture. He argued that the zeros of Ricmann's zeta function could bo put in 
one-to-one correspondence with the zeros of the s-wave Jost function for cutoffed 
potentials in the complex momenta plane. Joffily considered the s-wave Jost func- 
tion /+ (7) in the dimensionless argument 7 = kR {R is the cutoff of the potential) 
and the transformation 

z = -i—^ (18) 
by which the lower half of complex 7-plane is mapped onto the critical zeta line. 



7. The fractal potential of Wu and Sprung (1993) 
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Wu and Sprung obtained a local potential given by an Abel integral equation based 
on the smooth number of Riemann zeros below E, N{E) = ^ In -I- 1, considered 
as the scmiclassical WKB density of states below i? of a potential. Their potential 
y{x) = ^ scaled with respect to V'(O) = Vq (Wu and Sprung chose Vq = 3.100737r) 
is given by 



(19) 



As a function of the fitted Riemann zeros, y{x; N) shows fluctuations with respect 
to the smooth potential y{x) and therefore y{x;N) is a fractal curve. Using a box- 
counting- method analysis, Wu and Sprung obtained the result that the first 500 
Riemann nontrivial zeros can be fitted with the fractal potential y{x; 500) of fractal 
dimension 1.5. 

8. Okubo's two-dimensional potential (1998) 

Okubo introduced 

d'^ Q d i 

Hok{x, y; Pok) = + »/3oky^ + «(1 - l3ok)x— + - , (20) 

where /3ok is a real parameter. He gave arguments that this two-dimensional 

Lorcntz-invariant Hamiltonian may be relevant to the RH. Some cigcnfunctions 
of his two-dimensional Hamiltonian corresponding to infinite-dimensional represen- 
tation of the Lorentz group have many interesting properties. Especially, a relation- 
ship exists between the zero zeta function condition (ZZFC) that should be satisfied 
by any nontrivial zeta zero 



f 

Jo 



dr — = , (21) 

1 -I- exp r 
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for Re z G (0,1) and the absence of trivial representations in the wave function. 
The eigenfunctions of Hoi-cj) = X(j) with z = ^ + iX are of the form 

POO 

(j){x,y)= drr''''^ exp{ixT^~^°'^)g{T + yT'^'"') , (22) 
Jo 

where g is an arbitrary function which vanishes as its argument goes to oo. The 
choice 

50 (r + yr^-) = — ^ (23) 

1 + cxp(T + yTP°^) 

yields a function <j)Q which is an infinite-dimensional unitary reahzation of the 
Lorentz group S0(1,1). If the so-called ZZFC condition is applied to the repre- 
sentation space of S0(1,1) does not contain any singlet representation of the group. 

Berry and Keating commented that Okubo's hamiltonian is a more general, 
relativistic oscillator than their own quantized xp. 

9. Mussardo's potential (1997) ^'^ 

Mussardo's potential W is given by 

i-W 



(24) 



m=l 

where n is the Mobius function. 

Since Mussardo applies, like Wu and Sprung, the semiclassical method, his work is 
based on Abel's integral equation. All that we said with respect to the work of Wu 
and Sprung applies here as well. 

Mussardo also suggested a resonance experiment to implement the primality 
test up to a maximum number N depending on a cutoff eg- Sharp resonances 
in the transmission amplitude Tn{E) of plane waves impinging on the potential 
W{x;eo{N)) would indicate that the numbers N are prime. 

10. de Oliveira and Pellegrino (2001) 

These authors considered discrete Schrodinger operators of the type 

{Hvu)n = Un+1 + Un-1 + XV„U„ (25) 

defining Vo = 0,V^ = lifn-|-lisa prime number, and = if not. They have 
found a localization-delocalization transition as a function of the potential strength 
A. 



11. Meiny-body hamiltonians: Boos and Korepin (2001, 2002) 

Many-body hamiltonians are interesting more from the statistical physics point of 

view than the spectral one. Boos and Korepin, showed that particular correlation 
functions in quantum spin chains can be expressed in terms of the values of Riemann 
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zeta function at odd arguments (for even arguments, it is known that ({^n) B2n, 
the even BernoulU numbers). They put forward the following conjecture 

Arbitrary correlators of the XXX antiferromagnet are described as certain combi- 
nations o/ln2, the Riemann zeta function with odd arguments and rational coeffi- 
cients. 

An extension of this result to complex arguments could be interpreted as decay 
(decoherence) of the spin chain correlators. 

Wc recall that the XXX finite TV chain antiferromagnet is described since 1928 
by the following Heisenberg hamiltonian 

JV 

Hxxx = + ^^'f+i + - 1) ' (26) 

i=i 

where the tr's are Pauli matrices. 

12. Crehan (1995) 

Crchan used a classical theorem of Hardy and Littlcwood stating that the position 
of the nth prime on the critical line is limited by tn < oT^n [a = const) to get a 
theorem asserting that there is an infinite family of classically integrable nonlinear 
oscillators whose quantum spectrum is given by the imaginary part of the sequence 
of zeros on the critical line of the Riemann zeta function. In addition, Forrester 
and Odlyzko,^'^ related the distribution of zeta critical zeros to generalized Painleve 
fifth differential equation. 

13. Conclusion 

The spectral interpretation of the imaginary parts of the nontrivial Riemann zeros 
stimulated mathematical physicists to propose several quantum hamiltonians with 
spectra that could be useful in tackling with this century-old problem. Here the 
suggestions made in this area are gathered together in order to catch a global glimpse 
of these research facts. I provide a number of useful comments that hopefully will 
be taken into account in future possible experiments. The most exciting issue is the 
possibility of direct experimental access to negative phase shifts in the physics of 
trapped BECs. This cutting-edge research is of direct relevance for the link between 
prime numbers and black holes. I also hint that the rapidly developing PT quantum 
mechanics,^** a specific type of non-hermitic hamiltonians symmetric with respect to 
parity and time reversal discrete operations that display real spectra, will have an 
important contribution in this field. In fact, the susy hamiltonian pair of Pitkanen 
and Castro and collaborators are not hermitic but have real spectra. 
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